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Single-molecule experiments revealed that the replication rate of T7 DNA polymerase 
has a peak value at tension of about 5.5 pN. However, the mechanism leading to this peak 
has only been partially investigated. Here, we present a two-state model to investigate 
the effect of mechanical tension on DNA polymerase activity. The model consists of 
polymerase and exonuclease. The polymerase’s kinetic pathway has been simplified into 
Michaelis-Menten form with two steps. One is the conformational change from “open” 
to “close”. The other is the base-stacking. The results are in good agreement with 
experimental observations. We also predict that the tension on template is beneficial to 
the fidelity of DNA replication. 
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1. Introduction 


DNA polymerase (DNAp) motor is a protagonist in the propagation of genome. 
It works in a cyclic fashion and moves along a partly single-stranded (ss), partly 
double-stranded (ds) DNA chain. In each catalytic cycle, a single nucleotide will 
be added to the daughter strand. If the incoming nucleotide is complementary to 
that on the template strand, a base pair will form. Otherwise, it will be removed 
by exonuclease due to the fidelity of DNA replication. 13 

Crystal structure of T7 DNAp? has revealed that the polymerase resembles a 
half-open right hand. The “palm” subdomain forms a cleft flanked by the “finger” 
and “thumb” subdomains. The three subdomains hold the daughter-template DNA 
in S-shape and position the incoming dNTP for incorporation into DNA. The palm 
subdomain contains the catalytic site where the nucleotidyl transfer takes place. The 
fingers subdomain interacts with template and positions the incoming dNTP, while 
the thumb subdomain primarily binds the duplex DNA in a sequence-independent 
manner along the minor groove. 

T7 DNAp contains a 3’-5’ exonuclease domain which proofreads newly synthe- 
sized DNA and corrects mismatched base pairs to reduce the inherent low error 
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Fig. 1. The effect of mechanical tension on the net replication rate catalyzed by T7 DNAp. 
Circular dots and triangles represent our experimental data,!? the dotted line represents one 
previous model’s results,?° and the solid line is our model’s (rp). 


3,5-12 of 10-5 to 1077. On the rare occasion that a mismatch is incorporated 


ratio 
into DNA, the DNAp must unwind at least m = 8-9 base pairs and move the 
daughter strand from the polymerase site to the exonuclease site for editing.’ 

Single-molecule experiments! '4 have investigated the effect of mechanical ten- 
sion on T7 DNAp activity. The replication rate is strongly dependent on tension 
and has a peak at a critical tension of about 5.5 pN. The peak value approaches 
the rate of rate-limiting step in ensemble measurement,!° while the critical tension 
is located at the vicinity of the crossover force in the force-extension curves for 
dsDNA and ssDNA.!® 17 Previous theoretical models! 14: 18-21 could not quantita- 
tively explain the origin of the peak. For example, the dotted line? in Fig. 1 is one 
result of those models. We will show in this paper that this feature can be naturally 
explained by a two-state model, the prediction of which is shown by the solid line in 
Fig. 1. In our model, the exonuclease is considered as independent of tension, while 
the polymerase’s kinetic pathway15 22:23 will be simplified into Michaelis-Menten 
form with two steps. One is the conformational change from “open” to “close”. 
The other is the base-stacking. We will introduce Van der Waals-like potential to 
describe the conformational change before the formation of covalent bond between 
the “O” atom of the 3/-hydroxyl group of the daughter strand and the “P” atom 
of the a-phosphate of the incoming dNTP. 
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Fig. 2. The two-state model. State 1 represents the exonuclease, which is independent of tension. 
State 2 represents the polymerase, which has been simplified into a Michaelis-Menten form with 
two steps. kı represents the conformational change from “open” to “close”. kg is the base-stacking. 
Both steps are affected by the tension on the template. 


2. Two-State Model 


Figure 2 shows our two-state model. State 1 represents the exonuclease, while state 
2 represents the polymerase. We assume a mechanism exists in T7 DNAp: 


(1) In the case of correct incorporation, the potential of state 1 is AGT higher than 
that of state 2. 

(2) Once an incorrect incorporation occurs, the potential difference between two 
states will be exchanged, that is, the potential of state 2 is AG’ higher than 
that of state 1. It can be estimated that AGİ ~ mAG(0), where AG(0) is the 
free energy required to melt one base pair. 


In calculation of net replication rate, the mismatch situation is excluded. We 
assume the transition from state 2 to state 1 in single-molecule experiment!? is 
induced only by the force on the template. At steady state, the probability that 
DNAp is found in state 1 or 2 is given by 


1 
Pl = Ty emlAGO)—Aa(A/ Fat’ (1) 
and 
_ 1 : 
P2 = Tomo N-AGO| kat’ (2) 


respectively. Here kp is the Boltzmann constant, T the temperature, and Ag(f) = 
flass(f) — tas(f)], where ass(f) and xvas(f) are the end-to-end distances per base 
of ssDNA and dsDNA at tension f, and can be obtained from measurements of 
the stretching of DNA strands subjected to tension.!® 1!” Figure 3(a) shows the 
probabilities versus tension. 


2.1. State 1: Exonuclease 


According to the widely accepted editing mode,’ the melted daughter strand does 
not contact the template. The depolymerizing rate, thus, can be considered as 
independent of tension. Experiment data!? also shows that the exonuclease rate, 
Te, is saturated if tension is higher than 45 pN. We assume that re is independent 
of the tension and determined mainly by the rate of motor backward stepping. It 
can be estimated from experiment data? as re œ 30 bases/s. 
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Fig. 3. The probabilities that DNAp is found in (a) two states and (b) the predicted error ratio 
versus tension on template. 


2.2. State 2: Polymerase 


The pathway of DNA synthesis catalyzed by T7 DNAp has been revealed. 15 22 23 
Both DNA and the correct dNTP bind polymerase rapidly with high affinity, fol- 
lowing which the polymerase undergoes a slow, rate-limiting conformational change 
from “open” (Epe DNAn) to “close” (Eš e DNAn), in which the orientation of the 
nascent base pair will be adjusted optimally in preparation for the nucleophilic at- 
tack?4 between the “O” atom of the 3'-hydroxyl group of the daughter strand and 
the “P” atom of the a-phosphate of the incoming dNTP. The chemistry step is 
fast, and is followed by a reversed conformational change, which accompanies the 
base-stacking and leads to product release and translocation of DNA for the next 
catalytic cycle. 

For saturated DNAp and dNTP, as well as low PP; concentrations in single- 
molecule experiment, !? the kinetic DNA synthesis pathway for T7 DNAp15 22-23 
be simply described by the Michaelis-Menten form as Fig. 2, where kı represents 
the conformational change of palm from “open” to “close”, and kə is the base- 
stacking, which accompanies the reversed conformational change. The two processes 
are inevitably affected by the tension on the template. 


can 
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Fig. 4. The orientation of the nascent (“+1”) base pair. 0 is the angle between “+1” and “0” 
base pairs; D, the diameter of B-DNA; d, the distance between the “O” atom of 3’-hydroxyl group 
of the daughter strand and the “P” atom of a-phosphate of the incoming dNTP; and do is the 
sum of the Van der Waals radii of O and P atoms. 


2.2.1. Conformational change of palm from “open” to “close” 


As the result of evolution, T7 DNAp adjusts optimally the orientation of the nascent 
base pair in palm to form a covalent bond and achieve an optimal replication rate. 
In ensemble measurement, !” the DNA primer /template only consists of a synthetic 
25/36-mer. In single-molecule experiment,'? however, the measured chain contains 
10,416 base pairs and a constant force was applied on the template. The maximum 
replication rate!*: 4 implies that DNAp will work most efficiently at critical tension, 
fe, with £ss( fe) = £as( fe). We assume the orientation of nascent base pair (“+1”) 
in single-molecule experiment,!*? 0 shown in Fig. 4, will vary with the tension as 


[rss(f) = Las(f)], that is 
OCF) = Omax[Tss( f) = tas(f)|/das 3 (3) 


where das % 0.34 nm, is the distance between the adjacent sugar-phosphates of 
B-DNA. The distance between the “O” and the “P” atoms before the formation of 
covalent bond in Fig. 4 can be given by 


d(f) = 4/d? + [2D sin(8/2)]? , (4) 


where D ~ 2 nm and do ~ 0.342 nm, are the diameter of B-DNA and the sum of 
the Van der Waals radii of O and P atoms,?° respectively. We introduce the Van 
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to describe the conformational change of palm from “open” to “close”, that is, 
the potential of “open” state is —U(f) higher than that of “close” state. Because 
the dsDNA in palm is in an S-shape* whose distance per base approximates to 
that of ssDNA, the displacement at f direction in this step can be neglected. The 
conformational change from “open” to “close” in single-molecule experiment!® can 
be corresponded to that in ensemble measurement! by 


Sketa Bi (6) 
ki(f) = k-e T/T, (7) 


The electrostatic interactions in palm is considered as independent of tension and 
has been taken into account in €. 


2.2.2. Base-stacking 


The reversed conformational change from “close” to “open” 15 shown in Eq. (8) 
accompanies base-stacking, which will organize the adjacent sugar-phosphate units 
from S-DNA to B-DNA geometry, and has to do additional work Ag(f). 


k24 
E% èe DNAn Ep e DNAn41 - (8) 
2— 
If we assume that the additional work contributes to the activation energy required 
to reach a transition state from the “close” state, the rate coefficient for this step 
is 


kot (F) = kay (O)eA9)/Fa? , (9) 


where k24(0) is contributed mainly by the excess energy in base-stacking. The net 
rate of base-stacking in Fig. 2 is 


ko(f) = k2+ (f) — k2- . (10) 
The polymerase rate in Fig. 2 can be given by Michaelis Law: 
kı(f)k 


k1 + kil f) + ka(f) 


The average replication rate then can be given by 


(rp) = p2fp — Pire. (12) 


Mod. Phys. Lett. B 2007.21:1097-1105. Downloaded from www.worldscientific.com 
by NEW YORK UNIVERSITY on 02/01/15. For personal use only. 


The Effect of Mechanical Tension on DNA Polymerase Activity 1103 


Table 1. The fitted parameters with 
Eq. (12) to experiment data.1% 
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3. Discussion 


With the proposed parameters from bulk experiments? 15:22:23; m = 8, ki(fe) = 
300 s71, ko (0) = 1200 s7! and kp_ = 18 s~', we can use Eq. (12) to fit the 
experimental data.!3 The fitted parameters are listed in Table 1. The solid line in 
Fig. 1 shows that Eq. (12) is in good agreement with the experiment data.!° 

The assumed mechanism can also be used to explain error ratio. Once a mis- 
match is incorporated into DNA, DNAp must unwind m base pairs and move the 
daughter strand from the polymerase site to the exonuclease site for editing.? The 
probability that the 3’-hydroxyl group of the daughter strand still stays at the 
polymerase site for next cycle is the error ratio, which is about e7™1AG(0)+Ag(f)] 
at steady state. Figure 3(b) shows the predicted error ratio versus tension on tem- 
plate. The referenced m = 8 and fitted AG(0) = 2.1kpT are reasonable because the 
predicted error ratio for free DNA is consistent with the widely accepted data.!” 
The tension on the template, thus, is beneficial to the fidelity of DNA replication. 

Although we used a fitting parameter 0Qnax = 17°, Fig. 5(a) shows that the 
angle varies from —6° to 6° with increasing f. This range awaits testing. The 
replication rate at low tension, 130 bases/s, implies that the orientation to —6° 
may be constrained by the “close” palm. Figure 5(b) shows that the conformational 
change from “open” to “close” is a rate-limiting step for all tensions. Near fe, there 
is a plateau for kı at f € [4.5,6.5] pN, which may correspond to the situation in 
ensemble measurement.!° 

Unlike previous theoretical models,!* 14: 18-21 our model assumes that T7 DNAp 
does not work in the conformational change from “open” to “close”. This assump- 
tion is consistent with that where the dsDNA in T7 DNAp complex is bent into an 
S-shape by extensive interactions with the thumb and fingers.* In the reversed con- 
formational change, T7 DNAp has to do additional work because of base stacking. 


4. Conclusion 


Our model has for the first time explained the peak replication rate observed in 
experiment.'? The tension on the template not only affects the base-stacking, but 
also makes a great impact on the conformational change from “open” to “close”, 
which is described by a Van der Waals-like potential related to the distance between 
the “O” atom of 3’-hydroxyl group of the daughter strand and the “P” atom of 
a-phosphate of the incoming dNTP. The net replication rate of our framework is 
in good agreement with the experiment data!’ as shown in Fig. 1. If the tension 
on the template in single-molecule experiment!’ is between 4.5 and 6.5 pN, T7 
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Fig. 5. (a) @ is the orientation of the nascent base pair, Ag, the work to convert one ds residue 
into a ss residue, and U, Van der Waals-like potential, versus tension on the template. (b) Rate 
constants versus tension. kı represents the conformational change from “open” to “close”, and k2 
does base-stacking. fe is the critical value. Af ~% 2 pN, is the plateau of optimal replication. 


DNA polymerase works most efficiently and its rate approaches that of ensemble 
experiment. We can also predict that the tension on the template is beneficial to 
the fidelity of DNA replication. 
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